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Abstract —Optimal Power Flow (OPF) dispatches controllable 
generation at minimum cost subject to operational constraints on 
generation and transmission assets. The uncertainty and variabil¬ 
ity of intermittent renewable generation is challenging current 
deterministic OPF approaches. Recent formulations of OPF use 
chance constraints to limit the risk from renewable generation 
uncertainty, however, these new approaches typically assume the 
probability distributions which characterize the uncertainty and 
variability are known exactly. We formulate a Robust Chance 
Constrained (RCC) OPF that accounts for uncertainty in the 
parameters of these probability distributions by allowing them 
to be within an uncertainty set. The RCC OPF is solved using a 
cutting-plane algorithm that scales to large power systems. We 
demonstrate the RRC OPF on a modified model of the Bonneville 
Power Administration network, which includes 2209 buses and 
176 controllable generators. Deterministic, chance constrained 
(CC), and RCC OPF formulations are compared using several 
metrics including cost of generation, area control error, ramping 
of controllable generators, and occurrence of transmission line 
overloads as well as the respective computational performance. 

Index Terms —Wind power integration, power system eco¬ 
nomics, optimal power flow, wind power uncertainty, wind power 
variability, optimization methods, chance constrained optimiza¬ 
tion, distributionally robust optimization. 

I. Introduction 

A. Motivation 

The continually growing penetration of intermittent renew¬ 
able energy resources, e.g., wind and solar photovoltaic, is 
revealing a number of drawbacks in existing power system 
operational procedures that may limit the integration of these 
new resources. Wind generation is an intermittent and not fully 
dispatchable generation technology that imposes challenges to 
least-cost, risk-averse management of generation and trans¬ 
mission assets. One approach to address these challenges is 
strategic investments in more transmission and controllable 
generation capacity to enhance the system flexibility Q, 12. 
These investments are costly and subject to a variety of regu¬ 
latory and policy limitations. On the other hand, improving 
operating protocols may create additional flexibility in the 
existing system by replacing ad hoc deterministic policies 
for limiting system risk from intermittent generation with 
probabilistic formulations that account for intermittency in a 
principled manner. 

Manuscript received April 23, 2015; revised August 7, 2015. 

Miles Lubin is with the Operation Research Center, Massachusetts Institute 
of Technology, Cambridge, MA, 02142 USA (e-mail: mlubin@mit.edu i. 

Yury Dvorkin is with the Department of Electrical Engineering, University 
of Washington, Seattle, WA, 98195 USA (e-mail: dvorkin@uw.edu I. 

Scott Backhaus is with the Center for Nonlinear Studies, Los 
Alamos National Laboratory, Los Alamos, NM, 87545 USA (e-mail: 
backhaus@lanl.gov i. 


Historically, these deterministic policies were designed to 
account for less challenging deviations of load from its fore¬ 
cast value. They have performed very well for their design 
conditions where fluctuations in load are a small fraction 
of the total load, however, they are not expected to cost 
effectively manage risk when net-load fluctuations are large. 
Recent regulatory initiatives, such as Federal Energy Regu¬ 
latory Commission (FERC) Orders 764 ||3l and 890 H have 
identified the need for a new generation of operating protocols 
and decision-making tools for the successful integration of 
renewable generation. 

In this manuscript, we implement a distributionally robust 
chance constraint (RCC) optimal power flow (OPF) model and 
compare it with the deterministic OPF and chance constrained 
(CC) OPF models. The deterministic OPF is a typical short¬ 
term decision-making tool used by a number of utilities and 
its implementation in this work aims to give a reasonable 
benchmark for comparison. The chance constrained (CC) OPF 
limits the probability of violating transmission or generation 
constraints using a statistical model of wind deviations from 
forecast values, which are parametrized with a zero-mean 
Gaussian distribution 13, iSl. As compared to the CC OPF, 
the RCC OPF is a generalization that allows for uncertainty 
in the mean and variance of the wind forecast error. To 
demonstrate the effect of using probabilistic methods in an 
OPF, we compare these different OPF formulations in the 
setting of vertically-integrated grid operations, specifically, on 
a modification of the Bonneville Power Administration (BPA) 
system. 

B. Literature Review 

Traditional deterministic OPF models 121 dispatch control¬ 
lable generation using the central (most likely) wind forecast, 
i.e., they do not endogenously account for the variability and 
uncertainty of wind generation il. Exogenously calculated 
reserve margins and heuristic policies 13,1!3 are often used 
to enforce additional security requirements to ensure system 
reliability. However, these heuristic approaches are limited in 
their ability to produce cost-efficient solutions na. Cl. 

Recently, a number of transmission-constrained OPF and 
unit commitment (UC) models based on stochastic program¬ 
ming 0, Ha-Ha, interval programming Ha, na, chance 
constrained optimization ED and robust optimization na- 
II 20 II have been proposed for endogenous risk-averse decision 
making. A common drawback of stochastic programming 
and chance constrained optimization is a requirement for 
accurate statistical models of wind generation uncertainty and 
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variability. In practice, wind generation is modeled using the 
uni- or multivariate Gaussian distribution, which unavoidably 
results in solution inaccuracy ll22l . In addition, stochastic pro¬ 
gramming typically requires the generation of a relatively large 
number of scenarios leading to impractically large computing 
times, even for relatively large values of the duality gap E3i . 
If the number of scenarios is reduced by means of scenario 
reduction techniques ll24l . the monetary benefits attained with 
stochastic programming may reduce accordingly El, ESI. 

In contrast, interval programming and robust optimization 
models allow wind generation fluctuations within a certain 
range around a central forecast. By disregarding the likeli¬ 
hood of individual scenarios within this range, these methods 
result in an overly conservative solution as compared to 
stochastic programming ca. Although the conservatism of 
robust optimization can be reduced by adjusting the budget of 
uncertainty ll2^ or using dynamic uncertainty sets IZTl . there 
is no systematic methodology to choose the value of budget 
of uncertainty a priori. 

Here, we leverage recent work by Bienstock et al ll2T]l that 
formulated and implemented a CC OFF using a cutting-plane 
approach and demonstrated its scalability to large systems. 
Based on the CC OFF, Bienstock et al |I2T] also envisioned 
an RCC OFF that allowed for uncertainties in the mean and 
variance of wind forecast errors, but did not implement or 
test the RCC OFF. In this manuscript, we build upon the 
RCC OFF formulation from ED to reduce the inaccuracy 
of assuming Gaussian distributions for wind forecast errors. 
The RCC OFF is implemented in JuMFChance ll28l . an open 
source optimization package developed to model ordinary 
and distributionally robust chance constraints in OFF and 
other optimization settings. We use a modification of the 
BFA system to demonstrate the performance of the RCC 
OFF implementation on large-scale interconnections. Our case 
study shows that robustification of chance constraints can 
result in cost savings and is a more effective decision-making 
tool to mitigate real-time active power imbalances, extreme 
ramping of generators and transmission overloads. 

The rest of this manuscript is organized as follows. Sec¬ 
tion im reviews a deterministic OFF, restates (with a few 
modifications) the CC OFF of IJD, and demonstrates how 
robustifled chance constraints and are incorporated into the 
proposed RCC OFF model. Section [HI] describes the BFA 
test system and compares the performance of the deterministic 
OFF, CC OFF, and RCC OFF models on this test system. Fi¬ 
nally, Section |IV] presents conclusions and possible directions 
for future work. 

II. RCC OFF Formulation 

Section III-AI defines notations used throughout Section |II] 
In Section ITI-BI we review the CC OFF formulation from fTH 
and restate it with a few modifications. Section lTl-Cl discusses a 
model for distributionally robust chance constraints which are 
incorporated into an RCC OFF formulation in Section III-DI 

A. Notations 
B - set of buses 

r - index of reference bus (G B) 


C - set of lines 

Q - set of controllable generators 

Gb - subset (c Q) of generators located at bus b 

yV - subset of buses with wind farms 

Pmn - susceptance of line (m, n) 

fmn - real power flow over line (m, n), MW 

9b - phase angle at bus i 

Pi - output of controllable generator i, MW 

db - demand at bus b, MW 

pinin _ minimum output of generator i, MW 

p-max _ niaximum output of generator i, MW 

fmn^ - capacity of line (m, n), MW 

RUi - max ramp-up of generator i in the OFF period, MW/h 
RDi - max ramp-down of generator i in the OFF period, 
MW/h 

w[ - forecast output of wind farm at bus b, MW 
ujb{t) - actual deviation from forecast wl at time t, MW 
Oil - linear coefficient of cost for generator i, $/MW 
Ci 2 - quadratic coefficient of cost for generator i, $/MW^ 

5 - nonphysical auxiliary variable 
B - bus admittance matrix 

TTb - 5th row of the inverse of the admittance matrix, after 
excluding the row and column corresponding to the reference 
bus 

In the rest of this manuscript, bold symbols denote random 
variables. In particular, ujb models deviations ujb{t) within the 
OFF period, which drive random fluctuations in controllable 
generator injections Pj, bus phases 6b, and power line flows 
(described below). We denote the total deviation from the 
forecast as = X^bew ‘^b- In the CC OFF, the deviations iVb 
are assumed independent and normally distributed with zero 
mean and known variance In the RCC OFF formulation 
discussed in Section III-Dl this assumption is relaxed by intro¬ 
ducing symmetric intervals [—phtfib] and [cr^ — d-l,al + d1] 
for the mean and variance of bJb- 

B. CC OFF Formulation 

The CC OFF formulation is derived from the following 
single-stage deterministic OFF: 

min Y^(ci 2 P? + c^lPi) (1) 

p,9 ^ 

subject to 

BbnSn = '^p^ + w[ - db, V6 G B, (2) 

tiGB 

rnin ^ ^ ^ rnax \j- ^ n f'i\ 

P^ <P^<P^ , Vz € y, (3) 

fmn — ^mn{^m V^TTZ, 77-} G /I, (4) 

-/™r < fmn < v{m, n} G C, (5) 

where B is the \B\ x \B\ bus admittance matrix defined by: 

Pmm { 772 , 72 } G C 

12k:{k,n}eS 1 ^^"’ m = n ( 6 ) 

0, otherwise. 

In the deterministic OFF formulation in ([D-®, the control¬ 
lable generation set points pi are optimized to minimizing the 
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total cost of generation for the forecast wind generation 
and demand db subject to operating constraints on generators 
and transmission lines. In the presence of deviations Wf, from 
uj[, we model the proportional response of generators: 

Pt=Pi- Q;*0. (7) 

Here, > 0 is the participation factor for controllable 
generator i. When = 1, the response rule © guarantees 

that generation and load remain balanced, but does not limit 
the magnitude of the response of the generators or the resulting 
flow on the power lines. 

As shown in IITTI . the deterministic OPF in © can be refor¬ 
mulated as a CC OPF by introducing probabilistic constraints 
on and and modelling the participation factors ai as 
decision variables. We present the CC OPF formulation as 
follows: 

min ^ {caipl + var(r2)a^) -f cupi) (8) 

i&g 

subject to: 

Bbndn = '^ Pi + wl - db, ybeB (9) 

P™” <p^< pT"", Vt g a (10) 

fmn= Pmn(.0m- On), V{m, u} G £ (11) 

\fmn\<fZT, V{m,n}G£ (12) 

P{p, - Oa, < pr") > 1 - ez, Vz G 0 (13) 

P{p^ - > p™”) > 1 - ez, Vi G a (14) 

P{-na, < RU,) > 1 - £*, yiGO (15) 

P{flai < RDi) > 1 — Ci, G Q (16) 

P^fnin P ^m) P i^mn^ iPm 

V{m,n}G£ (17) 

P^fmn P ^m) P '^n) 

> > 1 - W, V{m, n} G £ (18) 

ai = \, a > 0 (19) 

i&G 

^ai=5r = 0r = Q. (20) 

ieGr 

Bbn5n = '^ cti, ybGB\{r} (21) 

uGB i^Gb 

n^r 

Here the decision variables are p, 9,6, a and /. For the 
quadratic cost of production (assumed to be convex), ED 
shows that 

IE[y^ Q 2 P- + caPi] = (qi(p- + var(r2)a^) -f caPi) ■ 
i&g ieG 

(22) 

Therefore, in ©, the CC OPF seeks to minimize the 
convex quadratic expected cost of production. Under these 
assumptions, the CC OPF is tractable and representable using 
second-order cone programming (SOCP). 


Constraints ©-(Ha are deterministic and enforce power 
flow feasibility, generator limits, and power line flow limits 
for u)^, similar to those in the deterministic OPF of ©. In 
chance constraints (fT3T l- (fT4l i. the controllable generator outputs 
are now random and given by Eq. ©. As in ED, these chance 
constraints bounded by Ci, the probability of the fluctuating 
generator outputs exceeding their upper limits p™“® or lower 
limits p™*". 

Relative to the CC OPF formulation of ED, we add only the 
constraints (flST l and (fThl l that limit the probability of the real¬ 
time response aiCl to wind deviations from the forecast from 
exceeding RUi for positive changes and RDi for negative 
changes. Here, RUi and RDi are the continuous ramping 
constraints on the generators over the OPF time step. 

Chance constraints (fT3T l- (fT6l l on the controllable generator 
injections are expressed explicitly in terms of a and the 
random wind fluctuations u). Chance constraints ([I7li-([l8]l, 
which bound the line flows are more subtle. The 

flows also change with the fluctuating wind injections 
and the controllable generator response, however, the 
depend on the wind deviations in an implicit manner, i.e. 

fmn — jdmniOrn ^n), where 

YBbnd„= 'ibGB. (23) 

n^B i^Gb 

Bienstock et al. ED derive explicit equations for by 
observing that once a reference bus r G Bis chosen and 9^ and 
ar fixed to zero, the system of equations (|2^ is invertible and 
the adjusted phase angles 9 (and hence /) can be expressed 
as a linear function of 9, fW, and a;: 

9b = 9b- flSb PTTbbJ (24) 

where 9 satisfies ®, S satisfies ED, and TTb is the feth 
row of B~^ (oriented as a column vector), where B is the 
(\B\ — 1 X |,B| — 1) submatrix of B with the row and column 
corresponding to the reference bus removed. The variable 6 is 
introduced solely for computational convenience. The chance 
constraints dnii-dn use EUl to express explicitly in 

terms of the random variables and decision variables. 

While chance constraints like (fT3l) - (fT8l) are often nonconvex 
and difficult to treat in general ll2^ . under the assump¬ 
tion of normality, they are both convex and computationally 
tractable EOl- In particular, a chance constraint of the form 

HGx <b)>l-e (25) 

is equivalent to 

p^x + $“^(1 — e)Vx'^T,x < b, (26) 

when 4 ^ £) where is the inverse cumulative 

distribution function of the standard normal distribution. In 
the following, we assume e < 4 so that — e) > 0 and 

constraint (l26l l is convex in {x, b). Note that in this model, we 
treat each chance constraint independently. Although it would 
also be natural to pose a model which attempts to enforce that 
multiple linear constraints hold jointly with high probability, 
convexity in this case remains an open question, even under 
the assumption of normality El, ED- 
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Constraint ( l26l l is not only convex; it can be represented as a 
second-order cone (SOC) constraint handled by many off-the- 
shelf optimization packages like CPLEX and Gurobi 1^ . 
Indeed, we see that (l26l l is satisfied iff 3t such that 

f>||S^x||2, (27) 

$"^(l-e)f < 6, (28) 

where (l27T i is a standard SOC constraint. 

While the representation of the CC OPF in (l8Tl- (l20l i with 
the reformulation of the chance constraints per (l26T l is quite 
useful, ED observed that off-the-shelf solvers were not ca¬ 
pable of solving large scale CC OPF instances. Instead, ED 
implemented a specialized algorithm based on sequential outer 
approximation of dZTl l. We note that in our implementation, the 
formulation is provided as stated in (I8ll- (l20l l to the modeling 
tool JuMPChance Il28l , which enables the user to select 
between solution via sequential outer approximation or via 
reformulation to SOCP. 

C. Cutting-plane algorithm to solve distrbutionally robust 
chance constraints 

In the analytic reformulation of the chance constraints 
in (|2^ , the wind deviations tUf, are assumed to be nor¬ 
mally distributed with known (zero) mean and variance, i.e. 
tu ~ E). This approach is computationally tractable 

but has drawbacks. The assumption of normality is often an 
approximation, and even when valid, p, and E are typically 
estimated from data and not known exactly. Often, we can only 
say with confidence that (^, E) fall in some uncertainty set (7. 
By reformulating the chance constraints of Section ITl-BI to so- 
called distributionally robust chance constraints, the constraint 

X < 6 is required to hold with high probability under all 
possible distributions within U, i.e., 

V(/r,E)GC/. (29) 

For each (/i, E) G U, we have a single convex constraint 
of the form (|25]) . therefore ( l29b is a potentially infinite set of 
convex constraints and is convex itself. In this manuscript, as 
in ED, we consider uncertainty sets U that can be partitioned 
into a product C7 = (7^ x Uy, where (/r, E) G C7 iff 
p G [7^ and E G Uy- Under these assumptions, Bienstock et 
al. ED observe that there is no known compact, deterministic 
reformulation of (|29] | and instead describe a cutting-plane 
algorithm which we build upon and demonstrate is capable 
of handling large-scale instances. In the rest of this section, 
our discussion departs from ED, highlighting a number of 
enhancements over the algorithm they propose. 

The cutting-plane approach in ED iteratively solves a se¬ 
quence of relaxations of E9] l. At each iteration, we must verify 
if EH* is satisfied. In the case of partitioned uncertainty, ED 
holds iff 

urwix^p, -f — e), /max x'^T,x < b. (30) 

For fixed x*, both of the inner maximization problems in ( |30] | 
have a linear objective, and so detecting if EH) is satisfied for 
a given x* can be computed by optimizing a linear function 



Fig. 1. An illustration of the iterative cutting-plane approach used to enforce 
the convex, nonlinear distributionally robust chance constraints. 

over the sets 77^ and Uy- If the solution to (l30l l shows that E9b 
is satisfied, then the algorithm terminates. Otherwise, dSOl l is 
used to find the corresponding (p*, E*) that violates ( l29l l, and 
we add a linearization of the corresponding constraint d26] l of 
the form 

x'^p* -\- $“^(1 — e)\J (a:*)'^E*a:*-|- 

- e)/^J{x*YE*x*'^ {x*fi:*{x -x*)<b (31) 

to the relaxation which cuts off the current solution. Therefore, 
at any iteration, the relaxation we solve is a linear program, 
similar to the linearization scheme for CC OPF. This process 
repeats until (1291 ) is satisfied within numerical tolerances. Our 
approach differs in two notable aspects from the algorithm 
proposed by ED: 

1) They propose to treat the term max^gj/ by stan¬ 
dard reformulation techniques based on strong duality 
to generate an equivalent deterministic formulation in an 
extended set of variables. Instead, based on the empirical 
observation that only a few extreme cases are important, 
we also apply the cutting-plane technique to this term. In 
other cases, it may be advantegous to reformulate E3- 

2) They propose to introduce a slack variable t for 
i/maxsgf/j, x^Ex and add a linearization whenever the 
constraint t > Y^maxEgc/j, x'^Ex is violated. However, 
it is possible for this constraint to be violated when the 
original EH is not. Our implementation avoids unnec¬ 
essary iterations and only adds linearizations when EH 
is violated. 

Although this algorithm does not have polynomial con¬ 
vergence guarantees in general, it is an immensely powerful 
approach that also underlies standard techniques such as 
Benders decomposition. See Figure [D for an illustration of 
the algorithm. We refer readers to ED for further discussion. 

We also note an enhancement for this algorithm based 
on the special structure of the power output and ramping 
constraints (HH-dlSl). For example, consider the constraint 

P{p, - na, < pTn > 1 - £*■ (32) 

Recall that Jl = U the total deviation from the 

forecast, hence 17 is a scalar Gaussian random variable with 
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mean and variance e^Se where e is the vector of In our RCC OFF formulation, we replace constraints (fTJt - 

all ones. Hence ( |32] | is equivalent to (fTSl l with their distributionally robust counterparts. 


Pi - fJ-bj ai + ^ ^(1 - e) Y^Q;2gT5]g < p->nax^ (- 33 ^ 

Vbew / 

In the CC-OPF case where we assume zero means and at > 
0, this constraint simplifies further to 

p, + $-1(1 - e)aye^Ee < (34) 

which is a single linear constraint; no SOC constraints are 
needed, an important property which was not noted by II2TII . 
Furthermore, in the distributionally robust case, the worst-case 
realizations may be computed independently of ai, i.e.. 


^(1—e)aj /max e^Se < 

V SgC/s 

(35) 

The constraint dTSl l is again a single linear constraint. This ob¬ 
servation is useful for constraints (fT3T l- (fT6l l with the structure 
where a single decision variable multiplies a scalar Gaussian 
random variable, since we may avoid applying the cutting- 
plane algorithm to these constraints. This special structure does 
not hold for the line capacity constraints (fTTI i- lfTsT l. 

Finally, we note that for a small number of special cases, 
including when U cannot be partitioned, Ben-Tal et al ll^ de¬ 
rive explicit semidefinite programming representations of (|2^ . 
We leave an exploration of this case for future work. 


Pi-at 



D. Formulation for RCC OFF 

We adapt the discussion in Section III-CI to formulate an 
RCC OFF First, we adopt the assumptions in ETll that the 
fluctuations at the different wind sites are independent within 
an OFF time step, i.e. S = diag/cr^). For brevity, we denote 
the multivariate Gaussian distribution with diagonal covariance 
matrix as where is a vector. The RCC OFF 

models uncertainty in the uj distribution parameters at each 
bus b as intervals [—/if,,/2b] and [cr^ — -\- al] for /r 

and a, respectively. To represent the aggregate uncertainty, 
we follow ED and construct polyhedral uncertainty sets 



for the mean and variance, similar to those proposed by 
Bertsimas and Sim 1^ . 


Pui^N(n,.^){Pi - < pn > 1 - ez, 


Vi €Q,yfi € Ufj,,a^ G Ucr2 

(38) 

Puir^N{ii,a^)iPi - ^ai > pT'") > 1 - et, 


Vi G C?, /r G C/^, cr^ G U„2 

(39) 

Puj^N{fi,(T'^)(. ClcXi ^ RUf) 6j, 


Vi G C?,/r G C/^,cr^ G U„2 

(40) 

Puj^N(fi,a^)i^^i — PPi) 


Vi G C?,/r G C/^,cr^ G U„2 

(41) 


(TTm “ 7r„) < > 1 “ 

V{to, n} G £, /r G cr^ G U ^2 (42) 


{^fmn 

“ 7r„) > > 1 “ 

V{m, n} G £, /r G C/^, cr^ G (43) 

For the objective (l22l i. we simply take the nominal value 
for var(r2), leaving this term unchanged. The parameters F^ 
and Fcr are the uncertainty budgets used to adjust the level of 
conservatism of the resulting RCC OFF algorithm. The F’s 
may be interpreted as a bound on the proportion of wind 
farms that may take their worst case distribution. The least 
conservative limit F = 0 recovers the standard CC OFF model. 
The most conservative limit F = 1 ensures feasibility when 
each wind farm can take on its worst-case production distri¬ 
bution. Uncertainty budgets have a well-studied probabilistic 
interpretation for the case of robust linear constraints 1^ . In 
this work, however, we enforce robust feasibility of a second- 
order cone constraint with respect to perturbations in the 
coefficients. One could interpret our formulation as presuming 
some prior distribution on the parameters of the distribution 
of UJ and then requiring feasibility of the chance constraint 
under high probability in the space of distributions. However, 
few probabilistic guarantees are known for this case; see ll34l 
Ch. 10]. 

Our motivation for choosing this formulation, although it 
lacks a rigorous probabilistic interpretation in terms of prior 
distributions, is that it adds very little computational expense 
to the CC OFF problem. Indeed, for the uncertainty sets in 
(l36l l and (ITTT i. the optimization (l30l l that required to check 
feasibility of each robust chance constraint (l29T l reduces to 
a cheap sorting operation EH, ESI- At the same time, as 
we intend to demonstrate in the case study, the RCC OFF 
formulation delivers measurable benefits in practice. 

Note that, similar to the CC OFF, the constraints (l38l l- (l43T l 
are treated separately, i.e. when they are “robustified” with re¬ 
spect to the distribution parameters, the worst-case distribution 
may be different for each constraint. In the remainder of this 
manuscript, we set F = F^ = Fc and investigate the results 
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as a function of F. 

We note a connection to the robust OFF models discussed 
in in ifTSl - l^ . ||25]| . IZTl . The reformulation from ( l26b to 
(I25I 1 under the assumption of Gaussian uncertainty has a 1- 
1 equivalence with robust linear constraints with ellipsoidal 
uncertainty sets 1261. Hence, our CC OFF may be interpreted 
as a particular form of robust OFF, and under that interpreta¬ 
tion, the effect of the RCC OFF formulation is to introduce a 
level of uncertainty sets on the nominal values of the robust 
problem. 

III. Case Study 

We investigate the performance and benefits of an RCC 
OFF approach relative to CC OFF and deterministic OFF ap¬ 
proaches by implementing all three in the setting of vertically- 
integrated grid operations, specifically, on a modihcation of the 
Bonneville Fower Administration (BFA) systerrQ- This setting 
allows us to evaluate the effect of using probabilistic methods 
in an OFF. We use F as a parameter to study the impact of 
robust conservatism on the cost of generation and the statistics 
of Area Control Error (ACE), generator ramping, and power 
line loading. 

A. Test System and Data 

We use a modification of the BFA systeirH with 2866 
transmission lines and 2209 buses including 676 load buses, 
176 controllable generators, and 24 wind farms. The total 
installed capacity of controllable generators is 40.6 GW, 
composed of 133 hydro generators (28 GW), 41 gas-fired 
generators (9.6 GW), 1 coal-fired generator (1.2 GW), and 
1 nuclear generator (1.8 GW). The total installed capacity 
of 24 wind farms is 4.6 GW. The technical characteristics 
of the controllable generators and the network configuration 
are adapted from Fower World 1^ . Tie-line power flows to 
neighboring interconnections are modeled as loads at the ends 
of the tie lines. 

The case study spans the Winter season from December 
2012 to March 2013 so that the factors affecting wind, 
load, and hydro generation are relatively stationary. Hour- 
ahead, hour-resolution load forecasts and their actual 5-minute 
realizations are taken from ll40l . As provided, these data are 
aggregated across the entire BFA system. We disaggregate the 
load forecasts and realizations among 676 load buses based 
on their population density. For each wind farm location, 
archived NOAA forecast data SD are used to provide hour- 
ahead, hour-resolution wind speed forecasts. Actual, 5-minute 
resolution wind speed realizations at each wind farm location 
are taken from BFA historical data ffOl . Both the forecasts and 
realizations are converted to wind power using an equivalent 
wind turbine for each wind farm location ll42l . We use the 
methodology from to compute the intervals and 

[cr^ — a'^+d^] for the uncertainty sets in Eq. (l36l l and dJTl l. 

The methodology in uses a data-driven statistical analysis, 

^We refer interested readers to ED for illustrative case studies of the CC 
OPF model on the IEEE Realibility Test System. 

^For visual representations of the BPA system and some of its partitions, 
interested readers are refeiTed to ED 


which reveals that hourly average wind speed is proportional to 
the intra-hour wind speed variability. Next, this relationship is 
used to fit historical wind speed forecast errors to a generalized 
normal distribution. This fitting yields the best-ht ranges on the 
hourly-average wind speed and its standard deviation. After 
that, these ranges are converted from wind speed to wind 
power units by using a typical wind turbine power curve. 

In this study, we seek to realistically emulate the short-term 
operational planning of the BFA system operator; however, not 
all required data are available. To hll in the missing data, we 
make the following operational assumptions: 

a) Nuclear Generators: The single nuclear unit in the 
BFA system is assumed to be a “must run” unit, and its hourly 
power output pi is set to the 95% of its nameplate capacity. 
Its participation factor in real-time balancing is set to zero 

HU. 

b) Gas and Coal Generators: The power outputs pi and 
participation factors ai of all gas-fired and one coal-fired 
generators are decision variables. 

c) Hydro Generators: Dispatch decisions for hydro gen¬ 
erators often depend more on water flow considerations rather 
than on power system operations ll44l . Instead of being co¬ 
optimized with thermal and nuclear generators, hydro dispatch 
levels Pi are fixed in all OFF formulations considered here 
as exogenous parameters ll45l . System aggregated, hourly- 
resolution hydro generation is taken from historical BFA data 
ll40i and disaggregated to individual hydro generators based on 
their installed capacity. The assignment of participation factors 
is also affected by water flow conditions beyond the scope of 
this work. Lacking operational data, we set the participation 
factors ai of the hydro generators to a common value, which 
itself is a decision variable. 

We refer interested readers to ll46l for the input data and 
the code used in this case study. 

B. Evaluation Procedure 

The evaluation procedure includes two steps, which emulate 
the hour-ahead scheduling and real-time dispatch, respectively, 
and are organized as follows: 

1) Step 1: The RCC OFF, CC OFF (i.e. RCC OFF with 
F=0), and a deterministic OFF are solved for hour t using 
hour-ahead, hour-resolution wind power forecasts uj{ and 
load forecasts d{ at each bus b and the generator commitment 
decisions ut^i described in Section IIII-AI The result is an 
optimal hourly dispatch p^^ and hourly participation factors 

2) Step 2: Next, these optimal decisions p^ ^ and a^i are 
evaluated in a quasi-static power flow simulation of the system 
behavior using the actual, 5-minute realizations of wind power 
0Jt^b{T) and demand dt^b{T) where r refers to the twelve 
5-minute intervals of hour t. For every t, we compute the 
actual power output of each controllable generator pt,i{T) as 
Pt.i(T) = max[p™”,min[pt,i(T),p™“^]] where 

PtA^) = Pt,i - X! (44) 

few 

Using ptAr), a DC power flow calculation yields the 
ft,mAA- The actual power output ptAA *^hen used to 
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calculate the actual hourly operating cost Ct by summing the 
cost of over all 5-minute intervals r. This emulation 

process reflects the vertically-integrated utility setting, i.e. 
costs are allocated according to energy delivered to the system 
with no markup cost for providing regulation. 

Steps 1 and 2 are repeated for every operating hour in the 
period from December 2012 to March 2013. The results for 
each T are analyzed for power system area control error (ACE) 
statistics, generator ramping statistics, and power line flow 
statistics as described below. 

C. Cost Performance 

The RCC OFF includes three user-defined parameters re¬ 
lated to its probabilistic nature. The first of these is T which de¬ 
termines the budget of uncertainty in wind forecast probability 
distributions defined by the uncertainty sets in (l36l l and dJTl l. 
Wind conditions change frequently, and we expect that this 
uncertainty parameter will be determined by short-term policy 
decisions of the system operators. In contrast, the parameters 
Ei and €mn limit the probability that equipment constraints are 
violated, i.e. generation output or ramping limits described in 
(fT3l l- (fT6l l and power flow limits on lines described in (fTTI i- 
These parameters are directly related to the impact on 
power system assets, and we expect these are determined by 
long-term, e.g. seasonal policies. 

In our case study, these parameters are set in a sequential 
process with the results shown in Tables U and |II] Fixing 
the parameters Ci and Emn at 1/6 and 0.0025, respectively. 
Steps 1 and 2 from Section IIII-BI are executed to determine 
the actual cost of ^neration (7“ over the entire study period 
for r between O.Cjj and 1.0 and for a deterministic OFF. 
The participation factors are not decision variables in the 
deterministic OFF. Instead, a fixed participation factor a = 0.05 
is used for all thermal and hydro generators. TableUshows the 
lowest (7“ is found at r*=0.6. 

Fixing r = r* = 0.6, Table HIl displays the sensitivity of the 
actual cost (7“ of the RCC OFF solution to the parameters a 
and Emn- We begin the discussion with the results in the row 
for emn=0.01. For this larger Cmn, potential network conges¬ 
tion plays less of a role, i.e. the selection of the generators’ 
Pj j and a/ ^ is less dependent on their location in the network. 
Instead, their selection is more sensitive to generator costs and 
constraints. Therefore, as decreases (moving left to right) in 
the emn=0.Ql row, the main impact is to spread the and 
a/j more uniformly across the fleet to reduce the ramping 
duty of any particular generator. As Ci —1 /48, more duty is 
placed on higher cost generators driving up the operating cost 
via the var(r2) term in ([8]l. However for our case study, the 
high percentage of very flexible hydro generators suppresses 
this cost increase—a result that is not expected to carry over 
to other power systems with different controllable generation 
fleets. 

Next, we consider the small Cmn limit shown in the 
er7m=0.0001 row in Table HIl Here, avoiding potential network 
congestion plays a larger role in the selection of the p^ ^ and 
a/j with the results becoming less sensitive to the generator 

^Recall that the CC OFF is the RCC OFF with r=0. 


cost and risk parameter e^. This is reflected in the elevated 
and relatively flat cost even as —?> 1 /48. In between the two 
extremes of Cmn, there is a relatively strong trade off in (7“ 
between Cmn (network risk) and (generator risk). 

The results in TablesHlandHIlsuggest that, for the test system 
used in our case study, the RCC OFF model achieves the best 
cost performance with F* = 0.6, = 0.0025, e* = 1/6. In 

the remainder of this manuscript, we assume that the long-term 
policy parameters are fixed at = 0.0025 and e* = 1/6, 
and present the technical analysis for variable F. 


TABLE I 

Cost performance of the OFF models in the period from 
December 2012 to March 2013 



CC 

OFF 

RCC OFF 

OFF 

F 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

- 

C“, 

M$ 

112.2 

111.6 

110.9 

108.7 

112.7 

114.9 

115.6 

A, 

M$ 

- 

-0.589 

-1.296 

- 3.471 

0.459 

2.716 

3.401 

A, 

% 

- 

-0.524 

-1.155 

- 3.093 

0.409 

2.420 

3.032 


The actual realized generation cost C“ for the period from December 
2012 to March 2013 is computed in Step 2 using the dispatch es p* ^ 
and participation factors aj' ^ computed in Step 1 of Section Illl-BI 
The least cost solution is found for r=0.6 and is marked in bold. 
Also displayed are the changes in cost and fractional changes in cost 
relative to the r=0 case. In all of these cases, emn=0.0025 and 6^= 
and 1/6. 


TABLE II 

SENSITIVITY OF THE RCC OFF COST (%*) TO Ci AND Emn 




- T2 

- 24 

II 

^mn 

= 0.01 

1.379 

1.471 

1.563 

1.563 

^mn 

= 0.005 

0.827 

1.195 

1.379 

1.471 

^mn 

= 0.0025 

0 

1.379 

1.379 

1.471 

^mn 

= 0.001 

0.460 

1.011 

1.471 

1.471 

^mn 

= 0.0001 

2.849 

2.941 

2.941 

2.941 


Fercentage changes in actual generation cost relative to the 
ei = 1/6, emn = 0.0025 case. All cases use F = F* = 0.6. 


D. Technical Analysis 

1) ACE Performance: The Area Control Error (ACE) is 
computed for each 5-minute interval r in hour t as: 

ACEtir) = - wt^bir)) - ^pt,*(r). (45) 

bGB ieg 

Figure 12^ displays the cumulative distribution function (CDF) 
of all the ACEtij) in the study period, and Figure |2]i 
displays the average of ACEtij). Starting from the most 
conservative F = 1.0, the average of ACEtir) displays a slow 
but monotonic increase showing that F effectively controls 
the system’s technical performance. At F < 0.6, the average 
ACEtir) increases more rapidly to the CC OFF at F=0 
and displays a significant jump from the CC OFF to the 
deterministic OFF. By accounting for fluctuations in wind, the 
CC OFF outperforms the deterministic OFF in controlling the 
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ACE, and by progressively accounting for uncertainty in the 
parameters of the distribution describing the fluctuations, the 
RCC OPE outperforms the CC OPE 

In Eigure |2^), the difference between the CDEs for E = 
1.0, 0.8, and 0.6 is not very significant. As our measure of 
conservatism is relaxed further (i.e. F <0.6), the CDEs show 
a general increase in the frequency of ACE events of all 
sizes and the emergence of a longer tail of large ACEt{T) 
values. We also note that F ^ 0.6 is the value where the 
ACE statistics first begin to significantly deteriorate and where 
RCC OPE cost takes on its minimum value. Above F =0.6, 
little additional ACE control performance is gained for the 
additional cost. This analysis suggests that, in the setting 
of vertically-integrated grid operations, the RCC OPE with 
an appropriately chosen F will result in better compliance 
with the control performance standards (CPS) BtI at a lower 
operating cost. 

2) Ramping Performance: The RCC OPE also reduces 
generator ramp rate (RR) violations as compared to the CC 


Eirst, as the robustness of the RCC OPE increases, the number 
of generators affected by RR violations is reduced from 11 
(for the CC OPE, i.e. F=0), to 4 with F=0.2, and to 2 with 
F=1.0. Second, the number RR violations per generator is 
also greatly reduced as F increases. It is noteworthy that 
both effects can be observed even for a relatively small level 
of robustness, e.g. F = 0.2. As shown in Figure [25), the 
greatest number of 5-minute intervals when the robust chance 
constraints on ramping are violated, is observed with F = 0.2 
for generator 159. This translates into the fraction RiO.004 of 
the total number of intervals considered in this case study, 
which is less than e* = 1/6 enforced in the chance constraints. 
Combined with the results from Table H] this analysis suggests 
that the RCC OPE with an appropriately chosen F achieves 
lower operating cost and avoids indirect costs related to wear- 
and-tear effects on controllable generators. In contrast, RR 
violations for the deterministic OPE (not shown in Figure O 
are observed on 39 generators—24 of which experience RR 
violations in more than one 5-minute interval. 



Area control error (ACE(r)), MW 


b) Average ACE(t) for different OPF models 



Budget of uncertainty, V 


Fig. 2. a) Cumulative distribution function (CDF) of ACEtij) for different 
values of r=0.0 to 1.0. The CDF for the deterministic OPF is not shown 
for clarity of the Figure, b) Average of ACEt{T) for the same values of 
r as in a). Under the conditions of the study, the generators always have 
sufficient downward flexibility to avoid ACEt(r) < 0, i.e. overgeneration 
is not observed. On the other hand, undergeneration {ACEt{T) > 0) occurs 
for all models. All results are computed with = 0.0025 and e* = 1/6. 


Fig. 3. a) Histogram of the number of generator ramp rate (RR) violations 
per generator over the entire study period for the CC OPF (i.e. the RCC OPF 
with r=0). b) Same as a) but for the RCC OPF with F = 0.2, 0.6, and 1.0. 
All results are computed with = 0.0025 and e* = 1/6. 

3) Transmission Overload Performance: From the 
uJt^b{T), and dt^b{T), a power flow solution yields ft,mn{T) 
from which power line overloads are computed. Figure 01 
displays a histogram of the number of overloads per power 
line for the four most frequently overloaded lines in the RCC 
OPF for F = 0.0, 0.2, 0.6, and 1.0. Several other lines are 
overloaded during the study period, but these overloads only 
occur during one 5-minute period. Interestingly, F does not 
have a significant impact on the frequency of overloads for 
the lines in Figure 0] 

Among the most overloaded lines in Figure 0] the greatest 
number of violations of robust chance constraints on power 
flow limits is observed on line 1813 with F = 0. These 
violations are observed on the fraction f«0.013 of empirical 
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Line number 


Fig. 4. Comparison of overloads that are observed with the CC and RCC 
OFF models during more than one 5-minute interval. All results are computed 
with = 0.0025 and e* = 1/6. 



realizations, which is larger than eL = 0.0025 enforced in the 
chance constraint^ For other lines in Figure @1 the empirical 
number of violations are less than e*. = 0.0025. 

‘’J 

E. Computational Performance 

The computations were carried out with CPLEX 12.6 ll^ 
as an LP solver on a Intel Xenon 2.55 GHz processor with 
at least 32 GB RAM on the Hyak supercomputer system 
at the University of Washington m- All modeling was 
done using JuMPChance 1281 . a freely available extension 
for the JuMP lISTI modeling language. For CC OPF, we used 
the cutting-plane algorithm instead of the second-order cone 
reformulation. For RCC OPF, we used the algorithm described 
in Section III-CI As shown in Figure |5] the average wall- 
clock time for an RCC OPF instance was seconds with 
an increase to ~20 seconds for instances with F = 1.0. Such 
an increase is expected because of the larger number of cutting 
plane iterations required for constrained problems. 

This case study demonstrates that an instance can be solved 
within seconds on a system with 2209 buses, which is com¬ 
parable to real-life power systems. Therefore, the proposed 
formulation is likely to be compatible with requirements of 
existing commercially available short-term planning tools. 

E Remarks 

In the presented case study we made several assumptions 
which are specihc to the BPA test system. As a result of 

^The exact cause of violations on line 1813 cannot easily be explained by 
particular attributes of the model or the BPA test system, but the issue can 
be resolved via out-of-optimization corrections (45), mi used in practice. 


these assumptions, the proposed methodology may lead to 
different cost performance and technical results, if applied 
to other systems. However, we note that the BPA system 
takes advantage of highly flexible and low cost hydro gen¬ 
eration, which are capable of accommodating variability and 
uncertainty of wind power generation at a relatively low cost. 
In other power systems, which predominantly feature fossil- 
bred and thus are less bexible and less cost effective in 
accommodating buctuations of wind power generation, the 
proposed methodology is likely to be of greater value. We 
refer interested readers to m, m, which extensively discuss 
the impact of generation mix’s bexibility and production cost 
on power system operations. 

IV. Conclusions 

Based on ED, we have developed and implemented algo¬ 
rithms to compute a distributionally-robust chance constrained 
optimal power bow (RCC OPF) that accounts for uncertainty 
in the parameters of statistical models that describe the de¬ 
viations of wind (or other intermittent) generation from its 
forecast. 

We have demonstrated the scalability of the RCC OPF by 
performing a seasonal case study on a modibcation of the BPA 
system. In this setting of vertically integrated grid operations, 
the case study shows that, compared to deterministic or even 
chance constrained OPF (CC OPF), the RCC OPF disttibutes 
both generation and regulation in a manner that can result in 
both cost savings and better technical performance; including 
fewer violations of transmission line limits, generator ramping 
limits, and smaller Area Control Error values. 

The work in this manuscript points to several areas for 
potential future work, including: 

• similar case studies should be performed on power sys¬ 
tems that are dominated by fossil generation instead of 
hydro generation 

• the RCC OPE formulation should be extended to include 
the effects of reactive power on nodal voltage magnitudes, 
transmission line limits, and generator limits 

• to better model generator ramping consttaints, the RCC 
OPE should be modibed to a time-extended or look-ahead 
formulation consistent with the operation of modernized 
power systems 

> the time-extended robust chance constraints should be 
adapted to day-ahead planning tools, such as UC and 
security-constrained UC. Further algorithmic develop¬ 
ments may be needed to tractably solve such formula¬ 
tions. 

• the current formulation should be extended to market- 
based operations to incorporate the cost of procuring 
frequency regulation capacity. 
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